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Abstract 

Cyber-physical systems (CPSs) are man-made complex systems coupled with natural pro¬ 
cesses that, as a whole, should be described by distributed parameter systems (DPSs) in 
general forms. This paper presents three such general models for generalized DPSs that 
can be used to characterize complex CPSs. These three different types of fractional oper¬ 
ators based DPS models are: fractional Laplacian operator, fractional power of operator 
or fractional derivative. This research investigation is motivated by many fractional order 
models describing natural, physical, and anomalous phenomena, such as sub-diffusion pro¬ 
cess or super-diffusion process. The relationships among these three different operators are 
explored and explained. Several potential future research opportunities are then articulated 
followed by some conclusions and remarks. 
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1. Introduction 


It is well known that the CPSs with integrated computational and physical processes can 
be regarded as a new generation of control systems and can interact with humans through 
many new modalities [l|. The objective of CPS is to develop new science and engineering 
methods in which cyber and physical designs are compatible, synergistic, and integrated at 
all scales. Besides, as we all know, the DPSs can be used to well characterize those cyber¬ 
physical process |2|, l3| and the actions and measurements of the system studied are better 
described by utilizing the actuators and sensors, which was first introduced by El Jai and 
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Pritchard in and mainly focused on the locations, number and spatial distributions of the 
actuators and sensors. 

Moreover, in the past several decades, fractional calculus has shown great potential 
in science and engineering applications and some phenomena such as self-similarity, non¬ 
stationary, non-Gaussian process and short or long memory process are all closely related 
to fractional calculus P, 6, 7|. It is now widely believed that, using fractional calculus in 
modeling can better capture the complex dynamics of natural and man-made systems, and 
fractional order controls can offer better performance not achievable before using integer 
order controls [8|, [Oj . 

Motivated by the above arguments, in this paper, let G be an open bounded subset of 
R” with smooth boundary dQ and we consider the following fractional DPSs: 


Zt{x, t) + Az{x, t) = u{x, t) in G X [0, b], 


( 1 . 1 ) 


where 6 > 0 is a given constant, u is the control input depending on the number and the 
structure of actuators and A may be a fractional Laplacian operator, a fractional power of 
operator or a fractional derivative. 

The contribution of this present paper is to analyze the relationship among the fractional 
Laplacian operator, fractional power of operator and fractional derivative and try to explore 
the opportunities and research challenges related to the fractional order DPSs emerging at 
the same time. To the best of our knowledge, no result is available on this topic. We hope 
that the results here could provide some insights into the control theory of this held and be 
used in real-life applications. 

The rest of the paper is organized as follows. The relationship among fractional Laplacian 
operator, fractional power of operator and fractional derivative are explored in Sec. [H In 
Sec. [3l the emerging research opportunities of the fractional order DPSs with those three 
operators are discussed. Several conclusions and remarks of this paper are given in the last 
section. 


2. Three different types of operators 


In this section, we shall introduce some basic relationships among the fractional Laplacian 
operator, fractional power of operator and fractional derivative. For further information, we 
refer the readers to papers from to [36] in the reference section of the present paper and 
the references cited therein. 


2.1. Fractional Laplacian operator and fractional power of operator 

This subsection is devoted to the difference between the fractional Laplacian operator 
and fractional power of operator. For more details, please see [lO, 111, 1^, 13, 14, 15] and 
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their cited references. 

Let us denote (—A)“/^ the nonlocal operator (also called the fractional Laplacian oper¬ 
ator) dehned pointwise by the following Cauchy principal value integral 


(-A)“/V(x) 


C^P.V. 



fix) - f{y) 


dy, 


0 < a <2, 


( 2 . 1 ) 


where = 


2°r(l/2+a/2) 


' 0 f|r(-a/ 2 )| ^ constant dependent on the order a. Obviously, we see that the 

fractional Laplacian (—A)“^^ is a nonlocal operator which depends on the parameter a and 


recovers the usua’ 
operator, see 
result. 


La pla cian as a —?■ 2. For more information about the fractional Laplacian 
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19 


20| and the references cited therein. Now we have the following 


Theorem 2.1. f^J Suppose that (—A)"/^ is defined in L^(0,/) for a G (0,2). Then, the 
eigenvalues of the following spectral problem 


(_A)“/2^(a;) = A^(a;), x E (0, 1), 


where f E L^(0, 1) is extended to all H by 0 is 


An 


nvr (2 — a) 7 r 

1 41 



( 2 . 2 ) 


(2.3) 


satisfying 


0 < Ai < A 2 ^ ^ Aj ^ . 


Moreover, the corresponding eigenfunctions fn of An, after normalization, form a complete 
orthonormal basis in L^( 0 ,/). 


Note that the constant in the error term O (^) tends to zero when a approaches 2 and 
in the particular case, when a = 2, we see that An = {nn/l)'^ without the error term. 

However, for a positive operator A on bounded domain [0,/]. Suppose 0 < Ai < A 2 < 
• • ■ < An < • • • are the eigenvalues of A, {,^ 1 , ^ 2 , ■ ■ ■ ,^n, ■''} s-re the corresponding eigen¬ 
functions and fn (* = 1) 2, • • •) form an orthonormal basis of L^(0, /). Let (•, •) be the inner 
product of L^(0, 1). We dehne the fractional power of operator A as follows: 


A‘‘f{x) = Y,>-n(^n,f)Ux), feV(0,l), (2.4) 

n=l 


Then A^ (i = 1, 2, • • •) are the eigenvalues of A^. This implies that the two operators are 
different. 

Moreover, the work spaces of the two operators (fractional Laplacian operator and frac¬ 
tional power of operator A) are different. Before stating our main results, we hrst introduce 
two Banach spaces, which are specihed in HHr 
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16|. 
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For f2 C R” is a bounded domain, s G (0,1) and p G [l,cxo), we define the classical 
Sobolev space hF^’P(r2) as follows 13|: 


:= i / G LP{n) : ^ g X fl) i (2.5) 

i |a:-2/|p+* J 

endowed with the natural norm 

ll/ll»'...(n) |/(x)rrfx + ^ ^ (2-6) 

is an intermediary Banach space between U’iyi) and hF^’^(r2). When a non-integer s > 1, let 
s = m + a with m G N and a G (0,1). In this case, let f with |/3| = m be the distributed 
derivative of /, then the classical Sobolev space hF®’^(r2) defined by 


W*’P(fi) := {/ G W™’^’(fi) : Df^f G W‘"’^’(fi) for all f3 such that \(3\ = m} (2.7) 


with respect to the norm 

Il/Ik-P!) - (ll/IIE,».,.(n, + (2-8) 

is a Banach space. Clearly, if s = m is an integer, the space coincides with the 

Sobolev space W'^’^{Vl). 

Besides, let p{x) ~ with 5{x) = dist(a:, 12^). Define another space as follows: 

W;’^(D) := {/ G W*’^(D) : p{x)f{x) G L^'(D)} (2.9) 


with the norm 




\p{x)f{x)\^dx + 



\fix) 


n Jn 


X 


y n+sp ^ 1 


( 2 . 10 ) 

Actually, Wp’^^Q) is called nonlocal Sobolev space and we have Wp’^^Q) C hF^’P(r2) [l^, 1^ . 

By the Remark 2.1 in llj], for the fractional power of operator, we take the classical 
fractional Sobolev space as its work space. But for fractional Laplacian operator, we must 
take the nonlocal Sobolev space as its work space, which can be regarded as the weighted 
fractional Sobolev space. More precisely, for any element / G hF*’^(r2), since —A is a local 
operator and we do not know how f{x) approaches 0 when x —?■ dQ. Even if considering the 
space W^'^in);={f G : fU = 0}, we only know that the function / = 0 on the 

boundary and we do not know how / approaches 0. However, for the fractional Laplacian 
operator, it is a nonlocal operator and it is defined in the whole space. So, it provides 
information about how / approaches 0 as x —)■ dQ. In fact, from the definition of nonlocal 
Sobolev space, we know that —)■ 0 as x ^ 912, which dictates how / approaches 0 near 
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boundary. It coincides with the result of the Theorem 1.2 in 22|]. Thus, this is a signihcant 
difference between the fractional power of operator A = —A and the fractional Laplacian 
operator. 

Besides, it is well known that || A/||vi/».p(o) = ||/||vF“+ 2 .p(r 2 )- But for the fractional Lapla¬ 
cian operator (—A)", ||(—A)“/||i 4 /s.p(q) = ||/||w'>+2“.p(o) will not hold. By using Fourier 
transform, we h^e ||(-A)"/||vy;.p(n) = \\f\\^s+ 2 ^,p^^y 

Finally, by 1^, let X be a rotationally invariant stable process of index a G (0,1). Its 
symbol is given by ri{f) = —1/|“ for all / G R"". It is instructive to accept that rj is the 
symbol for a legitimate differential operator; then, using the usual correspondence fj —)■ —idj 
for 1 < j < n, we would write 


A = -di -ad =-(-Ar/=, 


( 2 , 11 ) 


In fact, it is very useful to interpret r^iVL) as a fractional power of the operator — A. However, 
for the fractional Laplacian operator, it can be defined as the generator of a-stable Levy 
processes. More precisely, if Xt is the isotropic a-stable Levy processes starting at zero and 
/ is a smooth function, then 




lim 

h^0+ 


E[f{x)-f{x + X,)] 
h 


( 2 . 12 ) 


This also indicates that the fractional power of operator —A and the fractional Laplacian 
operator are different. 


2.2. Fractional Laplacian operator and fractional derivative 

This subsection is focused on the relationship between the fractional Laplacian operator 
and the Riesz fractional derivative. 

Definition 2.1. fil/ The Riesz fractional operator for n — l<a<nona finite interval 
0 < X < I is defined as 

where 

1 Or /* 2^ 

= --rj— / (x - 

F(n — «) ox^ Jo 

and 

are the left-sided and right-sided Riemann-Liouville fractional derivative, respectively. 
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Moreover, according to 2^, the fractional Laplacian is the operator with symbol \x\^. In 
other words, the following formula holds: 

= (2.14) 

where and denote the Fourier transform and inverse Fourier transform of /(x), 
respectively. We refer the readers to |25[ for a detailed proof of the equivalence between the 
two dehnitions fl2.ip and fl2.14p of fractional Laplacian operator. 

By using Luchko’s theorem in 2^, we obtain the following result on the equivalent rela¬ 
tionship between the Riesz fractional derivative and the fractional Laplacian operator 
_(_A)“/2; 

Lemma 2.1. fj^/ For a function f{x) defined on the finite domain [0, /], and /(O) = /(/) = 
0, the following equality holds: 


(-aA/W = 


2cos(f^; 


[oDff{x)+^D^f{x)] = 


(2.15) 


where a G (1,2) and the space fractional derivative is a Riesz fractional derivative. 

For more information on the analytical solution of the generalized multi-term time and 
space fractional partial differential equations with Dirichlet nonhomogeneous boundary con¬ 
ditions, we refer the readers to j^. For more information on the numerical solution of 
ifferential equation with Riesz space fractional derivatives on a hnite do- 


fractional partial c 
main, consult 


28 


29|. 


2.3. Fractional derivative and fractional power of operator 

In this part, we hrst show the following dehnition of the positive operator. 

Definition 2.2. f!^/ The operator A is said to be positive if its spectrum (j{A) lies in the 
interior of the sector of angle (p G (0,7r), symmetric with respect to the real axis, and if on 
the edges of this sector, Si = : 0 < p < oo} and S 2 = {pe"®*^ : 0 < p < 00 }, and 

outside it the resolvent {XI — A)~^ is subject to the bound 

M(p) 


||(A/-kl)-^|| < 


(2.16) 


1 -h |A| 

Moreover, for a positive operator A, any a > 0, one can dehne the negative fractional power 
of operator A by the following formula 3l| : 


R-" = — / A-“R(A)dA, (R(A) = {A- XI)-\ F = U R 2 ) • 
2m Jr 


(2.17) 


It is then quite easy to see that A~°‘ is a bounded operator, which is an entire function of a, 
satisfying = A~'^ if a is an integer n, and = A~°‘A~^ for all a,/3 G C. Using 

(12.171) . we have 


= — / A"“R(A)dA + — 


2'Ki 


2 m 


A"“R(A)dA, 


(2.18) 


6 

























Then taking the integration along the lower and npper sides of the cnt respectively: A = se 
and A = it follows that 


nQ7r2 /*oo 


AI-" = 


2ni 


s °R{—s)ds 


-.—am roo 


27ri 


s °'R{—s)ds 


cos(a:7r) + isin(a7r) 


sin OTT 


TT 


2m 

noo 


s °‘R{—s)ds — 


coslavr ) — i smi^avr ) 
2'Ki 


s °‘R{—s)ds 


s °'R{—s)ds 

poo 

s~‘^R{—s)ds 


r(«)r(i - a) Jq 

Moreover, for any a G (n — 1,n), we get that 


A^f = A^-^A^f = 


r(n — Q;)r(l + a — n) 


s'^-^R{-s)A^fds. 


(2.19) 


and for more properties on the fractional power of a positive operator, please see 


32 


33 


M, 


35l | and the references cited therein. 

Now we are ready to state the following resnlts on the connection of fractional derivative 
and integral with fractional power of positive operator. 


Theorem 2.2. Let the absolutely space 

21C”|0,/| := {/ : 6 0(0,/|, /<”)(l) 6 L=l0,/|} 

and let A be the operator defined by the formula Af{x) = f'{x) with the domain {f : f E 
AC^[0fi], = 0}. Then A is a positive operator in the Banach space y4C"'[0,/] and 

A‘^f{x) = oDff{x), n — 1 < a < n (2.20) 


for all f{x) G D{A). 


Proof. By 3^ , the operator A + si (s > 0) has a bonnded inverse and the resolvent of A 
is given by 

((21 + s/)-‘/) (x) = r e-‘^’->\f{y)dy. (2.21) 

Jo 

Then the operator A is a positive operator in /] and Eq. (I2.19p gives 

1 


A-f{x) = 


r(n — Q;)r(l + 0 — 77.) Jo 
1 

r(n — Q;)r(l + 0 — 77.) Jo 
1 

r(77 — Q;)r(l + « — 77.) Jq 

_1_ r 

T{n — Q;)r(l + a — n) Jq 


s^-^R{-s)A^f{x)ds 

{A + sl)-^ f^^\x)ds 

px 

/ e-^^^-AfA)fy)dyds 
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Let s{x — y) = X. we get that 


Then 


/ = {x- / X^-^e-^dX = {x - |/)"-"-^r(a - n + 1), 

'0 io 


AV{x) = 


r(n — a)r(l + a — n) Jq 
1 




r(n-a) 7o 

= oD:f{x). 


X_y^n-a-lfin)^y)dy 


This completes the proof. 


3. The emerging research opportunities 

Recent advances in modeling and control of fractional diffusion systems, fractional reaction- 
diffusion systems and fractional reaction-diffusion-advection systems have been reviewed in 
the framework of CPSs. The fractional order DPSs have now been found wide applications 
for describing many physical phenomena, such as sub-diffusion or super-diffusion processes. 
At the same time, to our best knowledge, many problems are still open calling for research 
cooperation of multi-disciplines such as mathematical modelling, engineering applications, 
and information sciences. 


First, it is worth noting that in the more recent monograph 37|, the theory of pseudo¬ 
differential operators with singular symbols, and the connections between them and those 
three types of operators are explored. See |38|, l39|, for more knowledge on pseudo¬ 


differential operator. Moreover, we claim that those equivalences between fractional Lapla- 
cian operator and fractional derivative, fractional order of operator and fractional derivative 
discussed in this paper can introduce new mathematical vehicles to study fractional order 
generalized DPSs. For example, when we study a fractional DPSs with Riesz fractional 
derivative, by Lemma 12.11 the spectral representation methods can be used to characterize 
the solution of the dynamic system. Then we can study the existence of solutions, stability, 
controllability and observability of the system under consideration. 

Potential topics such as modeling the sub-diffusion or super-diffusion processes with 
consideration of the networked mobile actuators and mobile sensors, the communication 
among the actuators and sensors, collocated or non-collocated actuators and sensors, their 
robustness and optimality problems are all interesting and worthy much more efforts in 
future. Another interesting and important topic concerns the time-space fractional DPSs 
where the traditional hrst order derivative is replaced by a fractional order derivative with 
respect to the time t. 

















Furthermore, in the case of diffusion systems, it is worth mentionin. 


not all the states can be reached in the whole domain of interest 


5 that, in general. 


41 


42| and it would 


be more challenging in nature since the dynamics of the real-life control problem is always 
hybrid continuous and discrete. Due to the spatial-temporal sampling and discrete nature 
of decision and control, the notions of regional analysis should be introduced, i.e., we can 
consider the regional stability, regional controllabili^ regional observability etc. of the 
system nnder consideration. In addition, as stated in from an application point of view. 


some plain qnestions such as “How many actnators/sensors are sufficient and how to best 
conhgure them for a fractional DPSs control process?”, “Given the desirable zone shape, 
is it possible to control or contain the fractional diffusion process within the given zone?”, 
if not, “how to quantify the controllability/observability of the actuators/sensors” and etc. 
might be asked, which in fact raises some important theoretical challenges and open new 
opportunities for fnrther research. 


4. Conclusion 

This paper is concerned with the fractional order DPSs with three different operators: 
fractional Laplacian operator, fractional power of operator and fractional derivative. The 
relationship among the three operators and the emerging research opportunities are intro¬ 
duced. We hope that the results here could provided some insight into the control theory 
analysis of fractional order DPSs in particular and CPSs in general. The results presented 
here can also be extended to complex fractional order DPSs and various open questions are 
still pending. For instance, the problem of regional optimal control of fractional order DPSs 
with more complicated sensing and actuation conhgurations are of great interest. 
Acknowledgement. This work was supported by Chinese Universities Scientihc Fund 
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